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Abstract 

Recently, a new model of propagation of the light through the 
so-called weakly three-dimensional Cole-Cole nonlinear medium with 
short-range nonlocality has been proposed. In particular, it has been 
shown that in the geometrical optics limit, the model is integrable and 
it is governed by the dispersionless Veselov-Novikov (dVN) equation. 

Burgers-Hopf equation can be obtained as H-l-dimensional reduc- 
tion of dVN equation. We discuss its properties in the specific context 
of nonlinear geometrical optics. An illustrative explicit example is 
considered. 

PACS numbers: 02.30.Ik, 42.15.Dp 
Key words: Nonlinear Optics, Integrable Systems. 

Many recent studies concerning the dispersionless integrable systems have 
shown their relevance in a broad variety of fields in physics such as Laplacian 
growth, topological field theory, nonlinear optics [l]-[in!, as well as their deep 
connection with various fields in mathematics, in particular, with the theory 
of asymptotic expansions, conformal and quasiconformal mappings and the 
study of integrable deformations of complex algebraic curves [TT]-[T^. 

In the present paper we are interested in the application of the theory of 
dispersionless integrable systems in nonlinear geometrical optics. In particu- 
lar, it has been shown in [TT that the Maxwell equations describing the propa- 
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gation of the light through the Cole-Cole weakly three-dimensional nonlinear 
medium admit an integrable geometrical optics limit. 

In this limit the system is governed by the dispersionless Veselov-Novikov 
(dVN) hierarchy, which is amenable by the quasiclassical 9-dressing method 

A reduction method based on the symmetry constraints jl5j appears to 
be an efficient approach to calculate solutions of the dVN equation [TB] . 

In what follows we will focus on a 1+1- dimensional reduction when the 
refractive index does not depend on one coordinate. In this case dVN hier- 
archy is reduced to the so-called Burgers-Hopf (BH) hierarchy. For sake of 
simplicity we will focus on the first equations of the hierarchy, the so-called 
dVN and Burgers-Hopf equations, respectively. 

In this regard, it is quite interesting to interpret some of the properties 
of BH equation within the specific context of nonlinear geometrical optics. 
One of them concerns the existence of breaking wave solutions. They could 
be useful to model dielectrics which present a kind of "impurities" inducing 
abrupt variations of the refractive index. Indeed, it can be seen that, at 
so-called breaking points, the curvature of the light rays blows up just like it 
happens at interface between different media. 

Let us start with the Maxwell equations in a dielectric medium 



VAH- — = V-D = 

dt 

dB 

VAE + — = V-B = 



(1) 



dt 

along with the constitutive equations 

D = eE, B = /iH. (2) 
Looking for monochromatic solutions 

E{x,y,z,t)=Eo {x,y,z) e"^"* 



(3) 



li{x,y,z,t) = Ho {x,y,z)e 
one gets from (0121) the following well-known second order equations 



V^Eo + uj'^fxeEo + (V log /i) A ( V A Eq) + V (Eq • V log e) = 
V^Ho + u; V^Ho + (V log e) A (V A Hq) + V (Hq ■ V log /i) = 0. 



(4) 

Medium under study is characterized by the following set of properties PlITU]: 
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1. e and /i obeys the Cole- Cole dispersion law 



e 

So + 



'^'7°''"' 0<.<i (5) 

We highlight that the range of values of the exponent v plays a crucial 
role in the construction of an integrable high frequency limit. 

2. sq and yUo depend only on the coordinates x, y and z, while e and jl 
are assumed to be depending on the coordinates, the fields and the 
spatial derivatives of the fields. The latter can be explained in terms 
of a mechanism of short-range non-locality by means of an integral 
constitutive relation among the electric field E and the displacement 
vector D 

3. All quantities, in high frequency limit show slow dependence on the 
variable z formally given by 

oz 

where ^ is a "slow" variable defined hj z = uj^^. 

Moreover, any of such a quantity /(x, y, z) can be expanded in asymp- 
totic series on the parameter uj'^ 

fix, y, z) = /(x, y, + ^'''fiix, y, + ^~^''f2{x, y,^) + .... 



Representing 

Eo = Eoe^"^, Ho = HqC^^^ (6) 

and using the previous assumptions in one of equations (jH), in high frequency 
limit (a; oo) one gets in the orders a;^ and a;^~^^ the main non-trivial 
contributions [Hj \W\ 

Si + Si = 4u, (7) 
Ss^ = (pix,y,^,S^,Sy) . (8) 

Equation (|7j) is the standard eikonal equation in two dimensions, cp is a. cer- 
tain function and 4u = eo/^o- As discussed in the paper [TU] the equations ((Tj) 
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and (jHI) constitute an overdetermined system for the phase S. The com- 
patibihty condition imposes specific restrictions on the possible forms of the 
function and the admissible refractive indices u. If is a polynomial differ- 
ential of S, the first non-trivial case is given by the third degree polynomial, 
i.e. 

^ = ^Sl - ^S,S'^ + V,S, + V,Sy. (9) 
The compatibility condition gives 

Vix - V2y = -3ux (10) 

+ Vly = 3Uy, 

which is the dVN equation. In the paper ^H] it has been shown how sym- 
metry constraints allow us to construct its 1+1-dimensional reductions of 
hydrodynamic type. A number of explicit solutions for the complex dVN 
equation have been discussed in (TU]. 

Particular situation in which refractive index depend only on one coor- 
dinate on the plane x-y is of physical interest too. So, let us assume that 

Uy = 0. (11) 

As a consequence of the eikonal equation the phase function S must 
satisfy the condition Sy = c, with c = const, then 

S = cy + S{x,0 (12) 

where S does not depend on y. Eikonal equation (jZj) becomes 

Si = 4u- c^ (13) 

while, choosing V2 = 0, one gets 

Vi = -3m, 

and the dVN equation (fTO|l is reduced to the Burgers-Hopf (BH) equation 

+ Quux = 0. (14) 
Solutions of BH equation can be expressed in the following hodograph form JTj 

X - 6< + ip{u) = 0, (15) 
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Figure 1: Typical breaking solution of Burgcrs-Hopf equation. Beyond breaking point 
> C*) inunction u becomes multi- valued. 



where i/j is an arbitrary function of its argument. Moreover, once calculated 
u{x,^) by the (fT^ the phase function is given exphcitly 

S^cy±jV*^d.. (16) 

Shock structure of the Burgers-Hopf equation is connected to the existence 
of a value ^* {breaking point), for which 

oo. (17) 

In the theory of partial differential equations these solutions are called break- 
ing waves. It is a textbook exercise to calculate the breaking point ^Tj. For 
instance, choosing the initial datum u{xo,0) = (1 — tanha;o)/6 the breaking 
point is ^* = 1 (see figure H}. 

We remark that, as discussed in ^H], the singular sector of BH hierarchy 
induces a stratification of the affine space of independent variables which 
gives rise to the integrable deformations of hyperelliptic curves. 

The propagation of the light on the plane vr : ^ = (^o = const) is 
governed by the standard eikonal equation ((7j). Denoting by S the congruence 
on 71 normal to the family of curves S{x,y,^o) = A, parametrized by A, the 
curvature of 7 G S can be represented by the well known formula [121 

K = — 1/ ■ Vm, (18) 
2u 

where V = {dx,dy) and u = {vi,i'2) is the unit principal normal to 7 on 
the plane vr. Now, let us observe that the curve 7 coincides locally with the 
projection of the light ray on the plane vr. 



Figure 2: Density plot of the real part of the refractive index n = rii + in^- Dashed 
line delimits the region where there is an imaginary contribution to refractive index (we 
denoted t = ^). 




Fi gure 3: -a- Three-dimensional visualization of the imaginary part of refractive index, 
-b- Three-dimensional visualization of the imaginary part of S. It vanishes everywhere 
except inside the absorption region. 
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Figure 4: Wavefronts configuration on the plane {x, ^). 



In virtue of the condition , the formula (fTHj) assumes the form 




(19) 



Then, if Ux blows up we conclude that the curvature of the light ray blows up 
as well. This is the typical behavior which light rays exhibit on the interface 
between different materials. This type of solutions could be useful to describe 
situations where light rays propagate in the non-homogeneous medium and 
at certain point they cross some kind of impurity. These impurities with 
drastically different optical properties could be responsible of the abrupt 
change of direction. 

As well known in electrodynamics, complex-valued refractive indices also 
have a physical meaning since they can be used to describe absorption effects 
of radiation. In this specific case, complex values of the refractive index n are 
associated with a strong damping of the electromagnetic wave. We illustrate 
this fact by the following explicit example. 

Let us consider a solution obtained from the hodograph relation (jlSp 
setting i(j{u) = u^. Expliciting u from (fT3jl and using it in the expression ()16|) . 
one gets 



We note that the refractive index n = I^Ju (we do not consider negative 
refractive indices) is not real-valued on whole plane x-^. Inside regions in 
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which the refractive index takes a complex contribution even the phase func- 
tion S becomes complex- valued. Thus, writing S = Si + 182, if 5*2 > the 
electric field acquires a damping factor 

E = Eoe-'^^^ e^'^^i. (22) 

In this regions strong absorptions effects are occurring. 

Figure El shows the real part of refractive index n = l^fu. The dashed line 
marks the region where the refractive index acquires an imaginary contribu- 
tion. This has been visualized in the three-dimensional plot in the figure Ol-a. 
Just where refractive index becomes complex-valued, the phase S takes an 
imaginary contribution indicating a strong absorption (see figure Ol-b). 
Figure m shows the wave fronts configuration in the plane x-^. Shaded region 
represents the absorption region. For ^ — cxo we have an almost plane 
wave front. Towards absorption region, wave fronts deform itself following 
level lines of refractive index. Just before absorption wave fronts present a 
deviation of an almost right angle. 
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